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\ Abstract. Recent results from fixed-target dimuon production experiments at Fermi- 

£^ ■ lab are presented. Various QCD tests using dilepton production data are discussed. 

We emphasize that clear evidence for scaling violation in the Drell-Yan process remains 
to be established. Further theoretical and experimental work are needed to understand 
the polarization of the Drell-Yan pairs. We also discuss the nuclear medium effects for 
dilepton productions. In particular, we discuss the use of Drell-Yan data to deduce the 
energy-loss of partons traversing nuclear medium. 
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^ : INTRODUCTION 
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The experimental detection of high-mass lepton pairs produced in hadronic reac- 
tions has a long and rich history. The famous quarkonium states that revealed the 
existence of the charm and beauty quarks in the 1970s were discovered through their 
dilepton decay branches. They are superimposed on a continuum, which was antic- 
ipated theoretically in 1970 [1], and is now known as the Drell-Yan (DY) process. 
The DY process, electromagnetic quark- ant iquark annihilation, is closely related to 
the deeply inelastic lepton scattering (DIS). By 1980, DY production was already a 
source of information about antiquark structure of the nucleon. Additionally, DY 
production with beams of pions and kaons yielded the structure functions of these 
unstable particles for the first time. Also notable in the history of the DY process 
were the discoveries of the W ± and Z° particles in 1983, produced by a generalized 
(vector boson exchange) quark-antiquark annihilation mechanism. 

New experimental work has been carried out in recent years by few but prolific 
collaborations working in the fixed-target programs at the CERN SPS accelerator 
and at Fermilab. A series of fixed-target dimuon production experiments (E772, 
E789, E866) have been carried out at Fermilab in the last 10 years. Some of the 
highlights from these experiments, namely the observation of pronounced flavor 
asymmetry in the nucleon sea and the absence of antiquark enhancement in heavy 
nuclei, are discussed by Garvey at this Conference. In this paper, we will focus on 
other areas of dilepton physics studied in these experiments. To follow the main 
theme of this Conference, we first discuss the status of various QCD tests using 
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FIGURE 1. Proton-induced Drell-Yan production from experiments NA3 [2] (triangles) at 400 
GcV/c, E605 [3] (squares) at 800 GeV/c, and E772 [4] (circles) at 800 GcV/c. The lines are 
absolute (no arbitrary normalization factor) next-to- leading order calculations for p + d collisions 
at 800 GeV/c using the CTEQ4M structure functions [5]. 



dilepton productions. We will then discuss the nuclear medium effects for dilepton 
productions. In particular, the relevance of the Fermilab experiments on the issue 
of parton energy loss in nuclei will be presented. 

QCD TESTS IN DILEPTON PRODUCTION 

In the "Naive" Drell-Yan model, the differential cross section is given as 

•:> d 2 a 8na 2 x±x 2 v-^ 2r , N , s , . , Nn /1X 
~dMdxF = —7^T 2 X (D 

Here q a {x) are the quark or antiquark structure functions of the two colliding 
hadrons evaluated at momentum fractions X\ and x-i. The sum is over quark flavors. 
The Feynman-rr (xp) is equal to x\ — x 2 - 

Although the simple parton model originally proposed by Drell and Yan enjoyed 
considerable success in explaining many features of the early data, it was soon 
realized that QCD corrections to the parton model were required. Historically, two 
experimental features demanded theoretical improvement: first, the experimental 
cross section was about a factor of two larger than the parton-model value, and 
second, the distribution of dilepton transverse momenta extended to much larger 
values than are characteristic of the convolution of intrinsic parton momenta. We 
now discuss several consequences of QCD corrections to the DY observables. 




FIGURE 2. Comparison between the NLO calculation from Ref. [6] and the E772 data [4]. 



Absolute Cross Sections and pr Distribution 

The inclusion of the NLO diagrams for the DY process brings excellent agree- 
ments between the calculations and the data. As an example, Figure 1 shows the 
NA3 data [2] at 400 GeV, together with the E605 [3] and E772 [4] data at 800 GeV. 
The solid curve in Figure 1 corresponds to NLO calculation for 800 GeV p + d (y/s 
= 38.9 GeV) and it describes the NA3/E605/E772 data well. 

Berger et al. [6] recently compared their NLO calculations with the E772 data. 
As shown in Figure 2, the px distribution is well reproduced by the calculation. 
At pt > 2 GeV/c, the DY cross section is shown to be dominated by processes 
involving gluons. This suggests the interesting possibility of probing gluon density 
using large pr DY events. 



Scaling Violation 

The right-hand side of Eq. (1) is only a function of xi,x 2 and is independent of 
the beam energies. This scaling property no longer holds when QCD corrections 
to the DY are taken into account. 

While logarithmic scaling violation is well established in DIS experiments, it is 
not well confirmed in DY experiments at all. No evidence for scaling violation is 
seen. As discussed in a recent review [7], there are mainly two reasons for this. First, 
unlike the DIS, the DY cross section is a convolution of two structure functions. 
Scaling violation implies that the structure functions rise for x < 0.1 and drop for 
x > 0.1 as Q 2 increases. For proton- induced DY, one often involves a beam quark 
with x 1 > 0.1 and a target antiquark with x 2 < 0.1. Hence the effects of scaling 
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FIGURE 3. Comparison of DY cross section data with NLO calculations using MRST [10] 
structure functions. Note that r = X\x<i- The E772 [4], E605 [3], and NA3 [2] data points are 
shown as circles, squares, and triangles, respectively. The solid curve corresponds to fixed-target 
p+d collision at 800 GeV, while the dotted curve is for p+d collision at y/s = 500 GeV. 



violation are partially cancelled. Second, unlike the DIS, the DY experiment can 
only probe relatively large Q 2 , namely, Q 2 > 16 GeV 2 for a mass cut of 4 GeV. 
This makes it more difficult to observe the logarithmic variation of the structure 
functions in DY experiments. 

Possible indications of scaling violation in DY process have been reported in two 
pion-induced experiments, E326 [8] at Fermilab and NA10 [9] at CERN. E326 col- 
laboration compared their 225 GeV n~ + W DY cross sections against calculations 
with and without scaling violation. They observed better agreement when scaling 
violation is included. This analysis is subject to the uncertainties associated with 
the pion structure functions, as well as the nuclear effects of the W target. The 
NA10 collaboration measured n~ + W DY cross sections at three beam energies, 
namely, 140, 194, and 286 GeV. By checking the ratios of the cross sections at 
three different energies, NA10 largely avoids the uncertainty of the pion structure 
functions. However, the relatively small span in y/s, together with the complication 
of nuclear effects, make the NA10 result less than conclusive. 

RHIC provides an interesting opportunity for unambiguously establishing scaling 
violation in the DY process. Figure 3 shows the predictions for p + d at y/s = 
500 GeV. The scaling-violation accounts for a factor of two drop in the DY cross 
sections when y/s is increased from 38.9 GeV to 500 GeV. It appears quite feasible 
to establish scaling violation in DY with future dilepton production experiments 
at RHIC. 
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FIGURE 4. Drcll-Yan angular distribution from Fermilab E772 [12]: p + Cu collisions at 800 
GeV/c. The dimuons cover the mass region 11 < M M + M - < 17 GeV/c 2 with —0.3 < xf < 0.8 and 
Pt < 6 GeV/c. Mean values for p t , xf, and M are 1.4 GeV/c, 0.16, and 11.9 GeV/c 2 , respectively. 
The solid curve is a fit to the data with the form 1 + Xcos 2 9, where A is 0.96 ± .04 ± .06. 



In the parton model, the angular distribution of dileptons is characteristic of the 
decay of a transversely polarized virtual photon, 



where 9 is the polar angle of the lepton in the virtual photon rest frame and A = 1. 
Early experimental data from both pion and proton beams [11] were consistent 
with this form but had large statistical errors. 

Recently, E772 has performed a high-statistics study of the angular distribution 
for DY events [12] with masses above the T family of resonances. About 50,000 
events were recorded from 800 GeV/c p + Cu collisions, using a copper beam 
dump as the target. Figure 4 shows the acceptance-corrected angular distribution, 
integrated over the kinematic variables. Analyzed in the Collins-Soper reference 
frame [13], the data yield A = 0.96 ± 0.04 ± 0.06(systematic). 

Including higher-order QCD corrections to the DY process [14,15] results in the 
more complicated form of the angular distribution, 
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where is the azimuthal angle and A, /i, and v are angle-independent parameters. 
NLO calculations predict [16] small deviations from 1 + cos 2 9 (< 5%) for p t below 



3 GeV/c. The relevant scaling parameter for the magnitude of these deviations is 
Pt/Q, implying that NLO corrections become important when p t — Q. A relation, 
1 — A — 2is = 0, developed by Lam & Tung [17], is analogous to the Callan-Gross 
relation in DIS. Measurements with pion beams at CERN [18] and at Fermilab [19] 
have shown that the Lam- Tung relation is clearly violated at large pt- 

Pion-induced DY experiments have unexpectedly shown that transverse photon 
polarization changes to longitudinal (A ~ —1) at large xp [18-21]. The xp depen- 
dence of A is qualitatively consistent with a higher-twist model originally proposed 
by Berger & Brodsky [22,23]. However, the quantitative agreement is poor. The 
model's basis can be described as follows. As x p of the muon pair approaches unity, 
the Bjorken-x (momentum fraction) of the annihilating projectile parton must also 
be near unity. Thus, the whole pion contributes to the DY process. This can be 
treated with perturbation theory, with the result that the transverse polarization 
of the virtual photon becomes longitudinal. The angular distribution at large xp 
becomes 

^ oc (1 -x) 2 (l + \cos 2 6) + asin 2 6 1 (4) 

where a is oc p 2 /Q 2 . 

Eskola et al [24] have shown that an improved treatment of the effects of 
nonasymptotic kinematics greatly improves quantitative agreement with the A val- 
ues from the pion data. Brandenburg et al [25] have extended the higher twist 
model to specifically include pion bound-state effects. They predict values for A, fi 
and v that are in good agreement with the pion data at large xp. Unfortunately, 
the results are quite sensitive to the choice of the pion Fock state wave functions, 
which are not well constrained by experimental data. 



NUCLEAR MEDIUM EFFECTS OF DILEPTON 

PRODUCTION 

From a high-statistics measurement of dilepton production in 800 GeV proton- 
nucleus interaction, the target-mass dependence of DY, J/\l/, and T productions 
have been determined in E772 [26-28]. As shown in Figure 5, different nuclear 
dependences are observed for different dilepton processes. While the DY process 
shows almost no nuclear dependence, pronounced nuclear effects are seen for the 
production of heavy quarkonium states. E772 found that J/\I/ and have similar 
nuclear dependence. The nuclear dependences for T, T' and T" are less than that 
observed for the J/\& and Within statistics, the various T resonances also have 
very similar nuclear dependences. 

Although the integrated DY yields in E772 show little nuclear dependence, it is 
instructive to examine the DY nuclear dependences on various kinematic variables. 
Using the simple A a expression to fit the DY nuclear dependence, the values of a 
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FIGURE 5. Ratios of heavy-nucleus to deuterium integrated yields per nucleon for 800 GeV 
proton production of dimuons from the Drcll-Yan process and from decays of the J/ip, ip', T(1S), 
and T(2S+3S) states [7]. The short dash and long dash curves represent the approximate nuclear 
dependences for the bb and cc states, A 096 and A 92 , respectively. 



1 .04 
1 .02 
s 1.00 
0.98 
0.96 



1.04 
1 .02 
B 1.00 
0.98 
0.96 







• 


• . • • 



0.1 



0.2 



8 10 12 14 

Mass 



0.0 



0.3 0.6 



1.0 2.0 



3.0 



FIGURE 6. Nuclear dependence coefficient a for 800 GeV p+A Drell-Yan process versus various 
kinematic variables [26]. 



are shown in Figure 6 as a function of xp(x 2 ), M, xp, and p t . Several features are 
observed: 

1. A suppression of the DY yields from heavy nuclear targets is seen at small x 2 . 
This is consistent with the shadowing effect observed in DIS. In fact, E772 
provides the only experimental evidence for shadowing in hadronic reactions. 
The reach of small x 2 in E772 is limited by the mass cut (M > 4 GeV) and 
by the relatively small center-of-mass energy (recall that x\x 2 = M 2 /s). p-A 
collisions at RHIC clearly offer the exciting opportunity to extend the study 
of shadowing to much smaller x. 

2. ol{xf) shows an interesting trend, namely, it decreases as xf increases. It is 
tempting to attribute this behavior to initial-state energy-loss effect. However, 
there is a strong correlation between xp and x 2 (xp = xi — x 2 ), and it is 
essential to separate the xp energy-loss effect from the x 2 shadowing effect. 
Figure 7 shows a versus xp for two bins of x 2 , one in the shadowing region 
(x 2 < 0.075) and one outside of it (x 2 > 0.075). There is no discernible xf 
dependence for a once one stays outside of the shadowing region. Therefore, 
the apparent suppression at large xp in Figure 6 reflects the shadowing effect 
at small x 2 rather than the energy-loss effect. 

3. a(p t ) shows an enhancement at large p t . This is reminiscent of the Cronin 
Effect [29] where the broadening in p t distribution is attributed to multiple 
parton-nucleon scatterings. It is instructive to compare the p t broadening for 
DY process and quarkonium production. Figure 8 shows A(p^), the difference 
of mean p\ between p-A and p-D interactions, as a function of A for DY, J/\l/, 
and T(IS') productions at 800 GeV. The DY and T data are from E772 [30], 
while the J/\I/ results are from E789 [31], E771 [32], and preliminary E866 
analysis [33]. More details on this analysis will be presented elsewhere [30]. 
Figure 8 shows that (pf) is well described by the simple expression a+bA 1 / 3 . It 
also shows that thep t broadening for J/ty is very similar to T, but significantly 
larger (by a factor of 5) than the DY. A factor of 9/4 could be attributed to the 
color factor of the initial gluon in the quarkonium production versus the quark 
in the DY process. The remaining difference could come from the final-state 
multiple scattering effect which is absent in the DY process. 

Baier et al. [34] have recently derived a relationship between the partonic energy- 
loss due to gluon bremsstrahlung and the mean p\ broadening accumulated via 
multiple parton-nucleon scattering: 

-dE/dz = \a. A(p 2 t ). (5) 

This non-intuitive result states that the total energy loss is proportional to square 
of the path length traversed by the incident partons. From Figure 8 and Eq. 5, 
we deduce that the mean total energy loss, AE, for the p+W DY process is ~ 0.6 
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FIGURE 7. Nuclear dependence coefficient a for the Drell-Yan process [26] versus xp for (a) 
x 2 < 0.075, right scale, and (b) x 2 > 0.075, left scale. The thin solid lines show a = 1. The 
dashed line is a linear least-squares fit to the lower points. Also shown is the mean value of x\ 
for (b). 
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FIGURE 8. The change of mean p 2 t for nuclear target, A(p 2 t ) = (p 2 )(A) - (p 2 )(D), for 800 
GeV p+A Drell-Yan process and J/* and T(15) productions. The solid curves are best fits 
to the A-dependence of A{p 2 ). The J/* and T(15) productions have identical curves for the 
A-dependence fits. 



GeV. Such an energy-loss is too small to cause any discernible effect in the x F (or 
xi) nuclear dependence. As shown in Figure 7, the dashed curve corresponds to 
AE = 2.0 ±1.7 GeV (for p+W), and the E772 data are consistent with Eq. 5. A 
much more sensitive test for Eq. 5 could be done at RHIC, where the energy-loss 
effect is expected to be much enhanced in A-A collision [35]. 
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